Microscopic study of spin-orbit-induced Mott insulator in Ir oxides 
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Motivated by recent experiments of a novel 5d Mott insulator in Sr2ir04, we have studied the two- 
dimensional three-orbital Hubbard model with a spin-orbit coupling A. The variational Monte Carlo 
method is used to obtain the ground state phase diagram with varying a on-site Coulomb interaction 
U as well as A. It is found that the transition from a paramagnetic metal to an antiferromagnetic 
(AF) insulator occurs at a finite U = Umi, which is greatly reduced by a large A, characteristic of 
5d electrons, and leads to the "spin-orbit-induced" Mott insulator. It is also found that the Hund's 
coupling induces the anisotropic spin exchange and stabilizes the in-plane AF order. We have further 
studied the one-particle excitations using the variational cluster approximation, and revealed the 
internal electronic structure of this novel Mott insulator. These findings are in agreement with 
experimental observations on Sr2lr04, and qualitatively different from those of extensively studied 
3d Mott insulators. 

PACS numbers: 71.30.+h, 75.25.Dk, 71.20.-b 
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Transition metal oxides have been one of the most fas- 
cinating classes of materials in recent years [l[ . For the 
past decays, tremendous amount of efforts have been de- 
voted to explore the nature of 3c? transition metal ox- 
ides where various exotic states and phenomena have 
emerged such as high-T c cuprate superconductors, colos- 
sal magneto-resistant manganites, multiferroics, and var- 
ious magnetic orders. It has been established that these 
states and phenomena are caused by strong Coulomb 
interactions along with cooperative interactions of spin, 
charge, and orbital degrees of freedom, which are basi- 
cally separable in 3d electrons 0. 

Very recently, 5c? transition metal oxides have attracted 
much attention as a candidate of a novel Mott insulator. 
Because of the extended nature of 5g! orbital, Coulomb 
interactions are expected to be smaller for 5c? electrons 
(~ 1-3 eV) than for 3c? electrons (~ 5-7 eV) 0, whereas 
the spin-orbit coupling (SOC) A is estimated to be con- 
siderably larger in 5c? (~ 0.1-1 eV) than in 3d (~ 0.01-0.1 
eV). Therefore, in 5c? transition metal oxides, inherently 
spin and orbital degrees of freedom are strongly entan- 
gled. 

One of such examples is the layered 5g! transition metal 
oxide Sr2lr04. In Sr2ix04, ti g and e g orbitals are well 
separated by the large crystal field, and the lower ti g 
orbital is filled with five electrons, fa g ) ■ In spite of 
the large band width and small on-site Coulomb interac- 
tion U, Sr 2 Ir04 is an antiferromagnetic insulator with a 
weak ferromagnetic moment 0, Hf. Neutron diffraction 
patterns do not detect any superlattice structure that in- 
dicates charge order or charge density wave states @ . It 
is proposed that the strong SOC is responsible for the 
insulating mechanism 0. Indeed, the Ad counterpart of 
Si'2Rh04, which has a larger U and a smaller A than 
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Sr 2 Ir04, is metallic @. 

The proposed picture of this 
Mott insulator in Sr2lr04 is as 
splits the t2 g orbitals into J c ff = 
±1/2, twofold degenerate) and J e g = 3/2 states (J* s = 
±1/2, ±3/2, fourfold degenerate). Here J e g denotes the 
"effective" total angular momentum derived from the 
large SOC with the large crystal field 0. When the SOC 
is large enough, the lower J e g = 3/2 state is fully filled 
and the upper J e g = 1/2 state forms an active half- filled 
energy band. The band width of this half-filled band 
(Wcff ) is much smaller than the original one without the 
SOC (W) as shown in Fig. QJa), and thus even small U 
can lead the system into a Mott insulator. This picture 
of "J c ff = 1/2 Mott insulator" is supported by several 
experiments such as resonant x-ray scattering [9j, angle 
resolved photoemission spectroscopy , and optical con- 
ductivity [7j. This novel Mott insulator has a quite dif- 
ferent character than 3c? Mott insulators where the effect 
of the SOC is only perturbative and thus the spin and 
orbital are essentially separate objects. Although this 
picture is supported by first-principle calculations based 
on the density functional theory 0, [TTJ] , a study treat- 
ing many-body effects beyond the mean-filed level is still 
lacking. The main purpose of this paper is to understand 
the nature of Mott insulator induced in spin-orbital en- 
tangled 5c? systems. 

Here, we study the novel Mott transition and mag- 
netic order induced by the SOC for Sr 2 Ir04. The ground 
state properties of the three-orbital Hubbard model with 
a SOC term are investigated with the variational Monte 
Carlo (VMC) method. We show that the large SOC 
works cooperatively with U and leads the system into 
a novel Mott insulating state with an in-plane AF order. 
We also calculate the one-particle excitation spectrum us- 
ing the variational cluster approximation (VCA) [12| to 
discuss the difference between the "spin-orbit-induced" 
Mott insulator and the 3g! Mott insulator. 
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FIG. 1: (a) Schematic picture of the splitting of density of 
states by the SOC. Non-interacting tight-binding energy band 
(b) without the SOC, (c) with the SOC (X/t = 1.028), and 
(d) with the SOC and staggered AF field. See the text for 
the tight-binding parameters used. 



We take the three-orbital Hubbard model on a two- 
dimensional (2D) square lattice defined by the follow- 
ing Hamiltonian H = -f/kin + ^so + Hi, where -f/kin = 
J2kaa £ a{k)c\ a<j c ka a is the kinetic term, H S o = A J2i L ^ 
Si is the SOC term with a coupling constant A, and 

Hj = U ^ n ia fn ia i, 

i,ot 

+ ^2 [VniavTlip-a + (U 1 - J)n laa nip,j} 

(i) 

is the Coulomb interactions including intraorbital, 
interorbital, and spin-flip and pair-hopping interac- 
tions [13l| . Here, a indicates electronic spins and the in- 
dices a and /3 represent three t% g orbitals, yz (1), zx (2), 
and xy (3). 

The kinetic and SOC terms can be combined (Ho = 
^kin + Hqo) m the matrix form, 

Ha = J2 ( c jL> c Lx> c L- ct ) 

ha 

( ei(fc) iaX/2 -crA/2\ / c kl<T \ 

x -iaA/2 e 2 {k) iA/2 c k2<J , (2) 

\-a\/2 -iA/2 e 3 (fc) / \c k3 - a J 

from which it is apparent that the SOC mixes the up- and 
down-spin states, and new quasiparticles are obtained 
simply by diagonalizing Hq. The new quasiparticles are 
characterized by the pseudo-orbital a and pseudo-spin a 
with a creation (annihilation) operator a kaa (a ka<7 ). In 
the atomic limit (e\(k) = E2(k) = £ 3 (fc) — 0), sixfold 
degenerate states are split into twofold degenerate J e g = 
1/2 states and fourfold degenerate J e s = 3/2 states 0j- 



First, we construct the non-interacting tight-binding 
energy band. In Sr2lr04, and also in Si'2Rh04, there 
is a large hybridization between the xy and x 2 — y 2 
orbitals originated from the tilting of IrOe octahedra. 
Because of this hybridization, the xy orbital is pushed 
down below the Fermi energy [l(| ■ To take into account 
the effect of this hybridization, we introduce the next- 
nearest and third-nearest hopping integrals in the xy or- 
bital. The chemical potential ^63 is also introduced to 
take account of both the hybridization and tetragonal 
splitting [Tl| . The resulting energy dispersion of the xy 
orbital is £ 3 (fc) = — 2£i(cos k x +cos k y ) — 4i 2 cosk x cosk y — 
2i 3 (cos2/c a; + cos2ky) + /x 3 . On the contrary, the yz 
and zx orbitals have almost one-dimensional character, 
£i(fc) = — 2t^ cos k x — 2t4 cos k y and£2(&) = — 2^4 cos /c^ — 
2t§ cos k y (£4 ^> £5). Using this form, the band dispersion 
of Sr 2 Ir0 4 calculated by the LDA+SO (spin-orbit) Q is 
well reproduced by choosing the tight-binding parame- 
ters with (£i,i 2 ,t3,t4,*5,M3> A)=(0.36, 0.18, 0.09, 0.37, 
0.06, -0.36, 0.37) eV, as shown in Fig. ^c). In the 
following, we set ii = t as an energy unit. To study 
the Mott transition and the role of the SOC, we only 
change the value of A and fix the other tight-binding pa- 
rameters. This assumption is justified by the fact that 
the LDA+SO band dispersion for the Ad counterpart of 
metallic Sr2Rh04 14] is well reproduced by choosing 
X/t ~ 0.5 with the other parameters fixed. 

Next, we examine the effect of Coulomb interactions. 
For this purpose, the following trial wave function is con- 
sidered, 

|*)=P Jc P G |cD>. (3) 

|$) is the one-body part obtained by diagonalizing 
Hq with renormalized (variational) tight-binding pa- 
rameters, Ho(ti, /i 3 , X a p). Note that by the effect of 
Coulomb interaction and tetragonal splitting, the "ef- 
fective" coupling constant of the SOC has orbital de- 
pendence: A — > X a p. To consider magnetically ordered 
states, a term with the different magnetic order param- 
eter is added to H . Here, we study the out-of-plane 
AF order (along z axis) and in-plane AF order (along 
x axis), described by A z £ e lQ - r '(al a ^a iat - a- Q ; a ^i) 
and A x £ e lQ T ' (a\ a ^a. lai + a\ ai a la ^), respectively. Here 
Q = (tt, 7r) is an ordering vector. Note that the staggered 
field is applied to newly formed quasiparticles in the real 
space (cii acr , a\ aa ) and a represents the pseudo-spin, not 
the original spin. The matrix to be diagonalized at each 
k becomes 12x12 for the AF state (not shown). 

The operator P G = [1 - (1 - 9j ) | 7 ) ( 7 |J in 

Eq. ([3]) is a Gutzwiller factor extended for the three- 
orbital system [l5j. Here, i is a site index and 7 runs 
over possible electron configurations at each site. For 
the three-orbital system, there are 4 3 = 64 electron con- 
figurations, namely, |0) = |0 0), |1) = |0 0t>, 
|63) = ItltltD- We control the weight of each elec- 
tron configuration by varying g 7 from to 1. The set 
of {g 7 } is a variational parameter and optimized so as 
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FIG. 2: (a) Ground state phase diagram of the three-orbital 
Hubbard model in a 2D 10x10 square lattice with n = 5, 
U' /U = 0.7, and J/U — 0.15. PM denotes the paramagnetic 
metal and a;-AFI denotes the AF insulator with an in-plane 
(along x axis) magnetic moment. The solid line indicates 
the first-order phase boundary Uya- The expected locations 
of Sr2lr04 (Ir) and S^RhCU (Rh) are also indicated in the 
phase diagram, (b) J/U dependence of Umi for different A. 



to give the lowest energy. In this study, we classify the 
possible 64 patterns into 12 groups by the Coulomb in- 
teraction energy, and g 7 's in the same group are set to 
be the same. 



E 



The remaining term Pj c = exp 

Eq. ([3]) is a charge Jastrow factor that controls the long- 
range charge correlation. The long-range charge correla- 
tion is known to be important for describing Mott tran- 
sition 16]. In this study, we assume that depends 
only on the distances, Vy — v(\r.i — For instance, 

we consider up to 19th-neighbor correlation for a 10x10 
square lattice, and therefore the number of independent 
variational parameters is 19. 

The ground state energy and other physical quan- 
tities are calculated with the VMC method. The 
variational parameters mentioned above are simultane- 
ously optimized by using the stochastic reconfiguration 
method (l7j . This method makes it possible to optimize 
many parameters efficiently and stably. 

Fig. Ufa) shows the obtained ground state phase di- 
agram for a 2D 10x10 square lattice with the electron 
density n = 5, corresponding to (^g) 5 , and the Coulomb 
interactions U' /U = 0.7 and J/U = 0.15. The transition 
from the paramagnetic metal to the AF insulator occurs 
at U = Umi, which depends sensitively on the value of A. 
This metal-insulator transition (MIT) is found to be first- 
order, indicating that the nesting of the Fermi surface is 
not perfect and the Coulomb interaction is essential in 
driving the transition. The insulating mechanism is un- 
derstood as follows. When the staggered (AF) field is 
applied, the degeneracy along the edge of the AF Bril- 
louin zone (along M-X in Fig. [IJ is lifted and the highest 
band is split off by the AF gap [Fig. H^d)]. As can be 
seen in Figs.QJb) and [He), the SOC lifts upward the two 
branches of the energy bands from the rest, which makes 
it easier to fully open the AF gap once the Coulomb in- 
teractions are considered. Namely, the larger the SOC 



is, the easier the system becomes the AF insulator. In- 
deed, Fig. [5Ja) shows that Umi monotonically decreases 
with increasing A. Moreover, we found that the effective 
coupling constant \ a p increases with U, indicating that 
U and A work cooperatively for insulating. 

We expect that this insulating mechanism found above 
is applied for Sr2ir04. Because the Coulomb interaction 
Uj T is much smaller than the band width, the metallic 
state is naively expected for Sr2ix04. However, Sr2ix04 
is experimentally found insulating with a canted AF or- 
der [9j. This counter- intuitive observation can be ex- 
plained if the SOC in Sr 2 Ir04, Ai r , is large enough to 
reduce Umi smaller than U\ r . Indeed, Ai r is estimated as 
large as 0.4-0.5 eV, which is much larger than the val- 
ues for 3c? and Ad electron systems. We consider that 
Ui r > Umi is satisfied and the "spin-orbit-induced" Mott 
insulator is realized in Sr2lr04. On the other hand, in 
Sr2Rh04, we consider that U^ih < Umi is satisfied and 
thus the metallic state is realized. The expected loca- 
tions of Sr2ir04 and Sr 2 Rh04 in the phase diagram are 
indicated in Fig. [2Ja), where both of them are located 
near the MIT point. Note that the observed insulating 
gap of Sr 2 Ir04 is very small (~ 0.1 eV), suggesting that 
the system is in the vicinity of the MIT point. It is also 
interesting to notice that a small amount of substitution 
of Ir (~ 10%) for Rh in Sr 2 Rh0 4 leads to the MIT [lj|, 
indicating that the system is located not far from the 
MIT point. 

The J/U dependence of Umi is also investigated and 
the results are shown in Fig. [2jb). It is clearly seen in 
Fig. (Hb) that Umi monotonically increases with increas- 
ing J/U, indicating that the Hund's coupling is unfa- 
vorable for the spin-orbit-induced Mott insulator. This 
behavior is naturally understood since the Hund's cou- 
pling competes with the SOC and works destructively 
for the formation of the quasiparticles originated from 
the SOC. It is, however, found that the overall shape of 
the phase diagram does not change qualitatively with in- 
creasing J/U except for Umi shifting to a larger value. 
For J/U =0.15-0.25, we estimate U M i =1.2-1.6 eV for 
Sr 2 Ir0 4 and 1.6-2.4 eV for Sr 2 Rh0 4 . 

In the insulating region, we found that the in-plane 
AF order (x-AFI) is more stable than the out-of-plane 
AF order (z-AFI). If there is no Hund's coupling [J = 
in Eq. ([T])] , the rotational symmetry in pseudo-spin space 
is preserved and z-AFI and x-AFI are energetically de- 
generate. However, the introduction of Hund's coupling 
induces the anisotropy and the in-plane AF order is more 
favored than the out-of-plane AF order. This result is 
consistent with the study of an effective strong SOC spin 
model by Jackeli and Khaliullin [19]. The magnetic x- 
ray diffraction experiment Q also supports the in-plane 
magnetism in Sr2lr04. 

We have also estimated the local magnetic moment 
as large as 0.3-0.4/j,b for the parameters appropriate 
for Sr2lr04. This value is comparable to the results of 
magnetic susceptibility measurements [13] and much 
smaller than the atomic value of 1 /j,b in the strong-SOC 
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FIG. 3: (color online), (a) One-particle excitation spectrum 
for the Mott insulating state. Spectra projected onto (b) 

Jeff = 1/2, (C) J off = 3/2,J c 2 ff = ±1/2, and (d) Jeff = 

3/2, Jeff' = ±3/2 states as well as (e) density of states p(to) are 
also shown. The parameters used are C/=1.44 eV, 7/'=1.008 
eV, J=0.216 eV, and A=0.432 eV, and w = corresponds to 
the Fermi energy. 



limit. This large reduction is due to the large itinerancy 
of bd electrons. If we assume that the magnetic moment 
exactly follows the tilting of IrOg octahedra observed ex- 
perimentally (~ 11°) [H, the expected ferromagnetic mo- 
ment is 0.05-0.07/iBi which is comparable to the experi- 
mental estimation [4j, |5j, [20| . 

Finally, to explore the internal electronic structure in 
the insulating state, the one-particle excitation spectrum 
is calculated using the VCA 21 1. The results are shown 
in Fig. [3] for a set of parameters appropriate for Sr 2 Ir04. 
As seen in Fig. GUJa), the Fermi energy is located inside 
the gap, and thus the state is insulating. The validity of 
the physical picture of "J c g = 1/2 Mott insulator" can 
be examined by projecting the spectrum onto J e g = 1/2 
[Fig. EJb)] and J cff = 3/2 [Figs. E^c) and EJd)] states 
defined in the atomic limit. It is clearly seen in Fig. [3] 
that the upper (unoccupied) band is well described by 



Jeff = 1/2 state. On the other hand, the lower (occu- 
pied) band is a mixture of J c fr = 1/2 and J G ff = 3/2 
states and can not be simply compared with the atomic 
limit. For example, the bands with J e ff = 3/2 character 
are located closer to the Fermi energy for all momenta 
except for the vicinity of the AF Brillouin zone boundary 
(M-X) [Figs. G3R) and (d)]. In this sense, the low-lying 
one-particle excitations are slightly different from the pic- 
ture of "J c [f = 1/2 Mott insulator" where the J c ff = 1/2 
and J e ff = 3/2 states are well separated and the former 
forms the lowest-lying upper and lower Hubbard bands. 
This difference is attributed to the large itinerancy of 5d 
electrons. 

Fig. [3je) shows the total density of states calculated 
using the VCA. As indicated by the arrow in Fig. [3]Je), 
the optical excitation is expected mainly at around 1 eV 
originating from ~ —0.5 eV to ~ 0.5 eV with a broad 
width due to a large band width of the occupied states. 
Indeed, a moderate peak structure is observed at around 
1 eV in the optical conductivity measurement Q , which is 
in good agreement with our calculations. In addition, the 
experiment has observed a rather sharp peak at around 
0.5 eV. This sharp peak may be explained by the ex- 
citation from the lower J c ff = 1/2 state to the upper 
Jeff — 1/2 state along the AF Brillouin zone boundary 
(M-X), indicated by the arrows in Fig. GUb), where the 
almost parallel bands expect to induce a sharp excitation 
(~ 0.5 eV). 

In conclusion, we have shown that the SOC plays 
an important role in determining sensitively the critical 
value of Umi as well as the nature of the resulting Mott 
insulating phase in S^IrCU. This is because the 5d Mott 
insulator, not like for 3d systems, is well described by 
the novel quantum number J e g, which is due to the large 
SOC along with the large crystal field, a generic feature 
for the 5d transition metal oxides. Therefore, we expect 
that not only U but also A are crucial factors to be con- 
sidered in describing MIT for the 5d systems in general. 
Moreover, because of the orders of magnitude different A, 
other novel quantum phenomena, such as the anomalous 
metallic state, multipole order, and unconventional su- 
perconductivity that are not observed in the 3d systems, 
are expected to emerge in the 5d systems. 

The authors thank Y. Yanase, M. Taguchi, and A. 
Rusydi for useful discussions. A part of this work is 
supported by CREST (JST). The computation in this 
work has been done using the RIKEN Cluster of Clusters 
(RICC) facility and the facilities of the Supercomputer 
Center, Institute for Solid State Physics, University of 
Tokyo. 



[1] M. Imada et al, Rev. Mod. Phys. 70, 1039 (1998). 
[2] Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 
[3] M. Taguchi, private communications. 
[4] M. K. Crawford et al, Phys. Rev. B 49, 9198 (1994). 



[5] G. Cao et al, Phys. Rev. B 57, R11039 (1998). 

[6] Q. Huang et al, J. Solid State Chem. 112, 355 (1994). 

[7] B. J. Kim et al, Phys. Rev. Lett. 101, 076402 (2008). 

[8] F. Baumberger et al, Phys. Rev. Lett. 96, 246402 (2006) 



■5 



[9] B. J. Kim et ai, Science 323, 1329 (2009). 

[10] It is a quite different band structure from the supercon- 
ducting counterpart of Sr2Ru04, where the :ry-originated 
7 band greatly contributes the superconducting behavior. 

[11] H. Jin et ai, Phys. Rev. B 80, 075112 (2009). 

[12] M. Potthoff et ai, Phys. Rev. Lett. 91 206402 (2003). 

[13] U — U' + 2J is assumed. See, e.g., J. Kanamori, Prog. 
Theor. Phys. 30, 275 (1963). 

[14] G.-Q. Liu et al, Phys. Rev. Lett. 101, 026408 (2008). 

[15] J. Biinemann et al., Phys. Rev. B 57, 6896 (1998). 

[16] M. Capello et al., Phys. Rev. Lett. 94, 026406 (2005). 

[17] S. Sorella, Phys. Rev. B 64, 024512 (2001). 



[18] H. Yoshida, private communication. 

[19] G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 102, 
017205 (2009). 

[20] N. S. Kini et al, J. Phys.: Condens. Matter 18, 8205 
(2006). 

[21] A \/2 x \f2 cluster is used. However, the results discussed 
here do not change qualitatively for different sizes of clus- 
ters. We have confirmed that the Mott transition is prop- 
erly described by the VCA. The fact that the two differ- 
ent methods (VMC and VCA) used give the same results 
indicates the validity of the spin orbit-induced Mott in- 
sulator. 



